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ABSTRACT 

- ^  A  mathematical  model  has  been  developed  to  predict  the  downstream  distribution 

of  particles  in  a  fluid  jet.  The  particles  were  initially  assumed  to  be  evenly  dispersed  over 
a  cross-section  of  the  circular  jet,  near  the  source  of  the  jet,  and,  with  computer-aided 
methods,  the  co-ordinates  of  the  particles  at  a  specified  distance  downstream  were 
determined.  The  computer  program  was  also  built  in  such  a  manner  that  the  paths  of 
specific  particles,  from  their  initial  upstream  position  to  their  final  position  downstream, 
could  be  completely  followed.  Some  graphical  results  of  these  paths  are  presented  for 
various  particle  sizes.  ^ — 
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NOTATION 


bi/i  half  the  width  of  a  circular  jet  at  half  depth,  cm 

d  particle  diameter,  cm 

D  jet  diameter  where  particles  are  initially  positioned  in  flow  stream,  cm 

00 

k  kinematic  momentum  of  circular  jet,  k  =  2n/0u2ydy 

L  jet  radius  where  particles  are  initially  positioned  in  flow  stream,  cm 

t  time,  seconds 

u  local  fluid  velocity,  cm  s"‘ 

u,  axial  component  of  fluid  velocity,  cm  s'1 

u,  radial  component  of  fluid  velocity,  cm  s"1 

U  center-line  velocity  of  fluid  at  any  distance  downstream,  cm  s‘‘ 

U0  center-line  fluid  velocity  at  x0,  cm  s*' 

x0  axial  distance  from  point  source  of  jet  where  particles  are  initially 

positioned  in  flow  stream,  cm 

x  co-ordinate  of  particle  position  parallel  to  center-line  flow  direction; 

axial  co-ordinate,  cm 

y  co-ordinate  of  particle  position  perpendicular  to  center-line  flow 

direction;  radial  co-ordinate,  cm 

a  particle  density,  g  cm"3 

M  absolute  viscosity  of  fluid,  poise 

q  fluid  density,  g  cm'3 

t0  virtual  kinematic  viscosity  of  turbulent  flow,  poise 

The  following  are  dimensionless: 

Co  drag  coefficient  for  spheres 

K  particle  inertia  parameter 

Re  spherical  particle  slip  Reynolds  number 

Re0  spherical  particle  Reynolds  number  in  free  stream 
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NOTATION  (Cont’d) 


time  tU/L 

axial  component  of  fluid  velocity  u,/U0 
radial  component  of  fluid  velocity  u,/U0 
dx/dT,  parallel  component  of  particle  velocity 
dy/dT,  transverse  component  of  particle  velocity 
parameter  used  in  calculation  of  fluid  velocity 
parameter  formed  by  combining  Re0  and  K 
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INTRODUCTION 

Previous  studies  have  been  made  (Mellsen,  1978  and  1979)  involving  streams  of 
suspended  particulates  in  such  a  way  that  impaction  forces  of  airborne  particles  and 
collection  efficiencies  of  specific  samplers  could  be  determined.  The  flow  considered  in 
the  present  study  is  of  this  same  type;  that  is,  it  involves  a  fluid  flow  with  suspended 
particles  moving  along  with  the  flow.  As  a  result,  the  same  mathematical  equations  as 
were  previously  used  may  be  used  here  to  describe  the  motion  of  the  particles  in  the  fluid, 
hi  doing  so,  the  velocity,  as  well  as  the  position  of  each  particle  at  any  distance 
downstream  can  be  determined. 

The  specific  flow  which  is  to  be  studied  here  is  that  of  a  turbulent  circular  jet  with 
particles  suspended  across  a  cross-section  of  the  jet.  By  developing  a  mathematical 
model,  the  radial  and  axial  co-ordinates  of  the  particles  at  a  specific  downstream 
distance  will  be  determined,  thus  the  actual  particle  distribution  will  be  known.  Also,  the 
model  is  to  be  devised  in  such  a  fashion  that  the  paths  of  particular  particles  may  be 
completely  followed  from  their  initial  upstream  position  to  their  final  downstream 
position. 
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The  information  as  revealed  by  this  model  will  prove  to  be  of  some  value  from  the 
chemical  defence  point  of  view  in  that  it  is  of  interest  to  know  how  particulate  agents, 
i.e.,  irritants,  such  as  CS  and  CR  disperse  in  a  flow  field  such  as  that  described  by  a 
circular  jet.  The  solution  may  also  be  applied  to  situations  involving  fluid-particle 
mixtures  flowing  from  smoke  stacks  and/or  exhaust  pipes. 

DEFINITION  OF  THE  PROBLEM 

When  a  jet  leaves  a  small  circular  opening  it  causes  surrounding  fluid  to  be  drawn 
in  and  mixed  with  the  jet  itself  (Fig.  1).  This  results  in  an  increase  of  volume  of  flow  with 
distance  from  the  orifice  and  a  downstream  cross-section  of  the  flow  field  which  is  much 
larger  than  a  cross-section  nearer  the  source  of  the  jet. 

The  problem  to  be  considered  here  involves  putting  particles  into  a  circular  jet 
flow  somewhere  upstream  and  determining  their  positions  downstream.  We  would  like 
this  situation  to  simulate  a  flu  id -particle  flow  from  a  tube,  with  particles  initially 
distributed  across  the  mouth  of  the  tube,  however  the  mathematics  describing  the  flow 
field  assumes  the  jet  as  coming  from  a  point  source.  In  order  to  account  for  this,  the 
particles  are  put  in  the  flow,  at  the  velocity  of  the  flow,  a  short  distance  downstream 
from  this  issuing  point.  Their  paths  are  then  followed  from  this  initial  point  onwards. 

The  distance  from  the  source  where  the  particles  are  initially  distributed  uniformly 
across  the  cross-section  is  defined  by  how  much  of  the  jet  velocity  profile  is  to  be  studied. 
It  has  been  found  that  a  distance,  xc  of  5.896D  from  the  point  source  encompasses  the 
width  of  the  jet  at  half  depth  (Schlichting,  1968).  That  is,  for  this  x0  value  the  velocity  of 
the  particles  which  are  located  on  the  circumference  of  the  circular  cross-section  will  be 
onc-half  that  of  the  center-line  velocity  (Fig.  2). 

By  increasing  x0,  a  smaller  amount  of  the  jet  angle  is  covered  and  the  velocity  of 
the  particles  on  the  circumference  will  be  more  than  one- half  of  the  center-line  velocity. 
When  the  value  of  x0  is  less  than  that  of  5.896D  the  opposite  situation  arises;  a  larger 
amount  of  the  jet  angle  is  covered  and  particle  velocity  on  the  circumference  becomes 
^smaller. 
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Thus,  the  basic  problem  is  to  choose  an  appropriate  x0  according  to  the  portion  of 
the  jet  angle  to  be  studied,  and  thereby  establish  the  starting  point  of  the  particle  paths. 
Having  done  this,  the  paths  are  to  be  traced  to  a  specified  distance  downstream  where  the 
co-ordinates  of  each  particle  are  to  be  determined. 

PARTICLE  MOTION 


The  motion  of  an  individual  spherical  particle  has  been  found  (Batchelor,  1967) 
to  be  defined  by  the  following  ordinary  differential  equations: 

--  H  (Eq.  1) 

dT  24  K 


dvx  _  CpRe(ux  -  vx) 
dT  24  K 


(Eq.  2) 


where  Re  = 

Re0((u,  - 

\yy  +  (ux  -  v,)1]*'1 

(Eq.  3) 

K  = 

odJU0 

18pL 

particle  inertia  parameter 
(Stokes  number) 

(Eq.  4) 

Re0  = 

U0de 

M 

free  stream  Reynolds  number 

(Eq.  5) 

The  variables  are  defined  on  the  notation  page^  at  the  beginning  of  this  report  and 
the  basic  geometry  of  the  flow  system  is  illustrated  in  Figure  1 . 

Both  Equations  1  and  2  involve  the  term  CDRe.  This  presents  a  slight  problem 
because  CD,  the  drag  coefficient  of  the  sphere,  has  been  found  to  be  a  function  of  the 
Reynolds  number  Re.  For  example, 


for  a  situation  involving  Stokes  flow.  However,  since  both  CD  and  Re  contain  a  velocity 
term  it  is  more  convenient  to  use 

Re  - 

24 


UNCLASSIFIED 


UNCLASSIFIED 


when  considering  Stokes  flow.  This  is  because  the  velocity  squared  term  contained  in  the 
denominator  of  the  definition  of  Cc  is  cancelled  by  the  velocity  squared  term  in  Re1  and 
C0Re  may  then  be  found  by  dividing  CcRe*  by  Re. 

As  a  result,  the  general  expression  becomes  Re  *  f(C0Rea)  and  in  actual  fact  the 
drag  coefficient  is  available  in  the  form  of  explicit  empirical  equations.  These  equations, 
are  (Davies,  1945): 


Re  =  - 


Co  Re* 


2.3363  x  lO-tCoRe1)1  +  2.0154  x  lO^CcRe*)’  -  6.9105  x  lO'^CrRe*)4 


(Eq.  6) 


for  Re  <  4  or  CDReJ  <  140 


log, o Re  =  -  1.29536  +  9.86  x  lO-Oog.oCoRe*)  _  4.6677  x  lO’fiog.oCoRe1)1  + 

1.1235  x  lO-'ClogCoRe1)*  (Eq.  7) 

for  3  <  Re  <  104,  or  100  <  CDRea  <  4.5  x  10’ 

The  assumptions  that  must  be  made  when  using  any  of  the  equations  which  were 
introduced  in  this  section  include  — 

(a)  negligible  effects  of  gravity  and/or  electrostatic  forces, 

(b)  monodisperse  spherical  particles  with  diameter  very  small  as 
compared  to  the  diameter  of  the  initial  upstream  cross-sectional 
area,  and 

(c)  free  stream  flow  that  is  steady,  incompressible  and  irrotational. 

THE  TURBULENT  CIRCULAR  JET 

The  velocity  components  of  a  turbulent  circular  jet  have  been  defined 
(Schliehting,  1968)  as: 


u„  = 


_3 _ k _ 1 

8n  t0x  (1  +  '/in*)3 


U,  =  —  \/3  —  — 5 
4  V  n  x  (1 

where  *7=4-  \f^~  ~~~ 

4  V  n  £0  x 


n  -  '/4*i3 
(1  +  'Arj*)* 


(Eq.  8) 
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The  curve  which  results  from  these  equations  has  been  found  to  agree  quite  favorably 
with  the  first  theoretical  consideration  of  a  circular  jet  as  done  by  W.  Tollmien  and  the 
experimental  results  of  H.  Reichardt  (Figure  24.8,  Schlichting). 

In  applying  Equations  8  to  the  problem  considered  here  it  is  necessary  to  find 
numerical  values  for  particular  constants  such  as  the  empirical  constant,  v/F/co,  and/or 
the  constant  kinematic  momentum,  k.  In  order  to  do  this,  the  theory  as  presented 
(Schlichting)  must  be  employed  to  modify  the  equations  of  u,  and  u,  to  more  suitable 
forms.  These  modifications  are  now  to  be  introduced. 

Due  to  measurements  obtained  by  H.  Reichardt,  half  the  jet  width  at  half  the 
depth  of  the  velocity  distribution  at  any  point  x  is  given  by  bi/2  =  0.0848  x.  This  is 
actually  half  the  width  of  the  jet  at  the  point  where,  for  each  distribution,  the  minimum 
velocity  in  the  x  direction  is  one- half  of  the  maximum  velocity  which  occurs. 

According  to  Schlichting,  the  point  at  which  the  fluid  velocity  is  one-half  the 
maximum  fluid  velocity,  rj  =  1.286  and  b,,2  =  5.27x  e0/vF.  By  equating  these  two 
expressions  of  b1/2,  the  result  is 

—  -  0.0161 
N/T 

so  that  \/k/£o,  which  occurs  in  u,  and  rj,  is  1/0.0161. 

The  equations  we  now  have  are: 

u,  =  —  (— !— ]  - ! -  Vf 

8nx  \  0.0161  (1  +  'Arf1)2 


u '  .  J  JT± 

4  v  n  x  (I  +  '/Srj2)1 

,  -  l  JT(-l-)  x 

4  V  n  0.0161  I  x 


It  can  be  seen  that  for  any  co-ordinate  consisting  of  an  x  and  y  value,  rj  can  be  found  and 
the  only  unknown  quantity  left  to  be  determined  is  n/F  found  in  ux  and  u,. 
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It  has  been  derived  from  substitution  that 
Vk  =  1.59b,,2U 

or  VTT  =  1.59  (0.0848  x)U  =  0.!348xU 

where  U  is  the  center- line  velocity  at  any  distance  x.  Not  knowing  this  value,  it  appears 
that  it  would  be  more  convenient  to  express  W in  tei  ms  of  U0  —  the  velocity  of  the  fluid 
at  x0,  where  x0  is  the  distance  downstream  from  the  issuing  point  of  the  jet  that  is 
considered  the  starting  point  of  the  solid  particles. 

It  is  known  that  for  a  free  turbulent  circular  jet  the  center-line  velocity  U  is 
proportional  to  x"1  therefore 

U  =  — 

x 

and  U0  =  —  where  c  is  an  arbitrary  constant. 

Xo 

By  obtaining  an  expression  for  x0,  then,  knowing  U0,  the  constant  can  be  found  as  can 
the  appropriate  value  of  U  for  any  x. 

Assuming  the  situation  as  illustrated  in  Figure  2,  then  b,/2  =  0.5  D.  Also, 
b,,2  =  0.0848  x0  as  was  previously  introduced.  Thus 

Xo  =  -M2_  =  5.896D  (Eq.  9) 

0.0848 

and  c  =  5.896DU,, 

The  velocity  U  at  any  downstream  distance  x  may  now  be  expressed  in  terms  of 
U0.  D  and  x: 

v  _  5.896DUq 
x 

and  the  simplified  expression  for  vT  becomes 
n/IT  =  0.795  DUo  . 

By  substituting  this  into  u,  and  u„  the  resulting  fluid  velocity  components  are: 

— 1 —  - ! -  (0.795  DUo) 

0.0161  (1  +  l/*n2)2 
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n 


u  y 


3  -  (0.795  DUo) 

(1  +  V*^)1 


Thus,  for  any  constants  D  and  U0  (which  are  analogous  to  OD  and  UF  of  the  computer 
program  used  to  solve  the  problem),  and  for  any  x  and  y  co-ordinate,  rj  may  be 
calculated  as  well  as  can  u,  and  u, . 

Each  fluid  velocity  component  was  normalized  in  the  solution  to  the  center- line 
fluid  velocity  which  prevailed  at  x0.  This  was  done  so  that  the  initial  velocity  in  the  x 
direction  of  the  fluid  at  the  center,  and  thus  the  center  particle,  was  always  1 .0.  All  other 
velocities  could  easily  be  compared  to  this  value. 


Finally,  as  an  overall  conclusion  to  these  previous  substitutions  and  derivations, 
the  normalized  equations  which  were  actually  used  to  find  the  velocity  components  of  a 
turbulent  circular  jet  for  the  situation  of  this  problem  are: 


u. 


u. 


3 

8nx 


1 

0.0161 


- - -  (0.795  D) 

(1  +  Arj2)2 


=  J-  \/X  _L  —21— .. ^ rl2_  (0.795 D) 

Uo  4  V  n  x  (1  +  'Art1)1 


(Eq.  10) 


where,  as  before 


COMPUTER  METHODS  OF  SOLUTION 


The  computer  program  which  was  built  to  solve  the  problem  considered  here 
consists  of  various  subprograms  incorporated  into  the  main  program  itself.  Each  of 
these  subprograms  are  briefly  described  in  the  program  listing  and  the  whole  program  is 
documented  in  such  a  way  that  the  step-by-step  procedure  may  be  quite  easily  followed. 
Refer  to  the  Appendix  for  a  listing  of  the  program  and  a  sample  solution. 

A  total  of  fourteen  constants  are  to  be  supplied  by  the  user  and  these  are  read  in  as 
input  data.  Each  of  these  variables,  as  well  as  most  of  the  others  used  throughout  the 
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program  are  described  at  the  beginning  of  the  listing  so  that  they  may  be  easily  referred  to 
while  running  the  program. 

In  order  to  establish  the  initial  upsream  starting  point,  a  constant  AF  is  to  be  read 
in  and  used  in  the  statement  XZ  =  5.896*OD/RADZ*AF.  The  value  of  A F  provided 
by  the  user  depends  upon  the  angle  of  the  jet  velocity  profile  to  be  studied.  With 
AF  =  1.0  the  jet  encompasses  the  jet  width  at  half  the  depth  of  the  profile.  See  Figure  2 
for  an  illustration  of  the  situation  resulting  when  AF  =  1.0  and  refer  to  the  previous 
section  entitled  “Definition  of  the  Problem”  for  a  more  complete  description  of  this 
value  of  x0.  When  AF  is  less  than  1.0,  then  a  different  value  of  x0  results  and  a  wider 
portion  of  the  velocity  profile  is  covered;  when  AF  is  greater  than  1.0,  a  smaller  jet  angle 
than  that  at  half  depth  is  covered.  Thus  by  varying  the  input  value  of  AF  different 
portions  of  the  jet  velocity  profile  may  be  studied  according  to  the  user’s  preference. 

One  of  the  most  important  subroutines  which  was  built  into  the  solution  program 
is  that  of  SUBROUTINE  SBM22.  As  is  described  in  the  listing,  this  uses  the  fourth- 
order  Runge  Kutta  method  to  integrate  a  specified  number  of  ordinary  differential 
equations.  In  this  case,  the  ordinary  differential  equations  which  are  to  be  solved  by 
SUBROUTINE  SBM22  are 

dv,  dv,  dx  dy 

IFF  ’IFF  'Hr  ’~dF 

which  are  analogous  to  DG(1),  DG(2),  DG(3)  and  DG(4)  respectively  of  the  computer 
program.  By  integrating  these  equations,  the  velocity  and  position  of  the  particle  is 
obtained  (note  that  DG(1)  and  DG(2)  represent  Equations  1  and  2). 

The  values  of  Re  and  K  of  Equations  1  and  2  are  easily  found  for  each  new  step  by 
direct  substitution  of  previously  determined  values  into  Equations  3, 4  and  5.  However, 
(he  value  of  CDRE  is  found  by  (he  subprogram  FUNCTION  CDRE.  This  is  done  by 
calculating,  for  each  step,  the  i.  merical  solution  of  the  definitive  empirical  equations 
(Eqs.  6  and  7).  By  using  Newton’s  method  for  finding  the  zero  of  a  function,  these 
numerical  solutions  are  determined  and  the  value  of  CDRE  can  be  obtained. 

This  subprogram,  FUNCTION  CDRE,  presents  the  only  visible  error  of  the 
whole  program.  The  procedure  used  involves  an  iterative  process  which  is  terminated 
only  when  a  certain  condition  is  met,  namely,  that 

x»*'  7  x»  «;  £ 

x*  ♦  i 
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where  t  is  any  small  number.  When  this  condition  is  not  met  then  a  message  of  NO 
CONVERGENCE  is  printed  by  the  computer  and  the  value  of  CDRE  is  taken  to  be  what 
was  previously  stored  there.  Although  the  number  of  iterations  (100)  is  quite  high,  there 
are  some  situations  where  the  initial  estimate  doesn’t  converge  and  the  message  is  printed 
out.  In  these  cases,  the  true  value  of  CDRE  is,  in  fact,  not  what  it  is  taken  to  be  and  a 
small  amount  of  error  is  introduced  into  the  calculations  by  assuming  an  incorrect  value. 
However,  since  the  actual  computed  values  of  CDRE  change  very  little  from  one  time 
increment  to  the  next,  then  this  error  is  assumed  to  be  negligible  and  the  message  may  be, 
for  practical  purposes,  ignored. 

The  program  was  written  for  a  large  number  of  diverse  situations  however  one 
limitation  which  was  discovered  is  that  the  smallest  particle  diameter  which  could  be 
input  is  20  m-  At  a  diameter  of  10  #i,  an  overflow  problem  in  the  computer  calculations 
resulted  so  no  solution  could  be  obtained.  Other  limitations  of  the  program  have  not  yet 
been  discovered. 

As  a  small  addition  to  the  solution,  the  average  number  of  particles  per  unit  area 
(cm2)  at  the  initial  cross-section  is  computed,  as  it  also  is  at  the  cross-section 
downstream.  It  may  be  noted  that  the  upstream  concentration  is  in  fact  the  actual 
number  of  particles  per  square  centimeter  since  they  (the  particles)  are  uniformly 
distributed  over  the  cross-section.  However,  upon  reaching  the  downstream  distance, 
the  particles  will  no  longer  be  uniformly  dispersed  and  the  resulting  concentration 
represents  an  average  of  particles  per  square  centimeter. 

RESULTS  AND  DISCUSSION 

A  sample  solution  for  one  set  of  input  data  is  shown  after  the  program  listing  in 
i he  Appendix.  Many  more  calculations  were  made  using  this  program  to  produce  the 
graphical  results  shown  in  Figures  3  to  7. 

The  plots  of  Figures  3,  4  and  5  are  those  of  the  actual  particle  paths  for  particle 
diameters  of  100  p,  30  and  20  \i .  It  can  be  seen  that  the  downstream  distribution  does 
not  vary  much  with  particle  size,  however,  the  larger  particulates  tended  to  become 
distributed  a  bit  further  away  from  the  center  of  the  flow.  Even  though  the  difference  is 
small  it  indicates  that  the  larger  particles  were  less  inclined  to  be  influenced  by  the 
surrounding  fluid  which  is  drawn  in  as  the  jet  spans  out  downstream.  If  this  is  the  case, 
the  momentum  of  these  larger  particles  would  cause  them  to  flow  out  and  away  from  the 
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center  more  than  the  smaller  particles  and  as  a  result,  they  would  have  greater  y 
co-ordinates.  It  must  be  noted  that  these  graphs  were  plotted  for  a  value  of  AF  equal  to 
0.3217  which  was  choosen  so  that  most  of  the  velocity  profile  was  covered  (this  was 
determined  for  v,  of  the  outside  paiticle  equal  to  0.0  at  the  initial  position). 

There  are  no  experimental  results  available  to  compare  with  the  mathematical 
results  presented  in  Figures  3,  4  and  5.  Since  the  model  contains  many  assumptions, 
including  steady,  irrotational  flow  and  the  starting  location  of  the  particles,  the  results 
should  be  tested  by  experiment  to  establish  their  validity.  Some  errors,  however,  could 
possibly  have  been  introduced  into  the  solution  by  computer  round-off  but  these  are 
assumed  to  be  negligible. 

The  graphs  of  Figures  6  and  7  are  families  of  curves  showing  the  radial  co¬ 
ordinate  of  the  outside  particle  as  a  function  of  the  inertia  parameter  K.  (The  outside 
particles  are  those  which  are  located  on  the  circumference  of  any  circular  cross-sectional 
area  of  the  jet  which  engulfs  all  particles.)  These  curves  are  plotted  for  a  series  of  values 
of  4,  a  dimensionless  parameter  independent  of  particle  size  and  equal  to  ReJ/K  . 

Figure  6  shows  the  results  for  AF  =  1 .0  which  means  that  all  of  the  particles  are 
assumed  to  be  distributed  across  the  width  of  the  jet  at  half  the  depth  of  the  velocity 
profile.  At  this  point,  as  was  previously  mentioned,  the  initial  velocities  of  the  outside 
particles  (in  the  x  direction)  are  one- half  the  velocity  of  the  center  particle.  Figure  7 
shows  the  results  for  AF  =  0.5  which  covers  a  wider  portion  of  the  jet  angle  than  when 
AF  =  1.0.  At  this  angle,  the  velocities  of  the  outside  particles  are  found  to  be  about 
0. 14  times  that  of  the  center  particle  thus  confirming  that  a  wider  portion  of  the  profile  is 
indeed  covered. 

In  these  two  figures  each  curve  was  fitted  by  eye  through  several  data  points 
obtained  by  applying  the  computer  program  to  various  particle  diameters.  As  a  result, 
by  knowing  the  particles’  free  stream  Reynolds  number  Re0  and  the  inertia  parameter  K 
then  by  studying  the  appropriate  curve,  the  radial  co-ordinates  of  the  outside  particles  at 
a  specific  distance  downstream  may  be  estimated.  Also,  for  any  value  of  $  which  falls 
between  those  plotted  in  the  figures,  the  method  of  interpolation  provides  a  means  of 
approximating  the  downstream  radial  co-ordinate. 
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As  a  last  point  to  note,  the  ret'  Its  as  given  in  the  Appendix  indicate  that  although 
both  the  particles  and  the  fluid  have  the  same  initial  velocity,  the  velocity  of  the  fluid 
tends  to  decrease  faster  than  the  particles  as  they  are  observed  downstream.  This  can  be 
explained  by  considering  the  effect  of  the  fluid  on  the  particles  in  the  following  way: 

In  the  presence  of  a  fluid  which  is  continually  decreasing  in  velocity,  as  is 
the  case  here,  particles  moving  in  the  same  direction  as  the  fluid  will  have  a 
drag  force  acted  upon  them.  Although  the  particles  will  begin  to  slow 
down  very  quickly  they  will  still  have  a  greater  velocity  than  the  fluid 
because  of  the  inertia  effects  which  they  experience.  As  the  distance 
downstream  is  increased  however,  the  difference  between  particle  velocity 
and  fluid  velocity  is  narrowed  (refer  to  columns  of  VX'  and  UX'  of  the 
solution  given  in  the  Appendix).  This  is  due  to  the  fact  that  the  fluid 
velocity  is  changing  more  slowly  as  the  downstream  distance  is  increased. 

CONCLUSIONS 

In  some  situations  it  may  be  of  interest  to  know  how  particulates  disperse  in  a 
fluid-particle  circular  jet.  As  an  attempt  to  find  some  answers  a  mathematical  model  has 
been  developed  which  predicts  the  distribution  of  particles,  by  specifying  their  radial  and 
axial  co-ordinates,  at  any  distance  downstream.  This  model  also  enables  the  paths  of 
specific  particles  to  be  completely  followed  from  their  initial  upstream  position  to  their 
downstream  position. 

A  computer  program  was  written  to  provide  quick  and  easy  solutions  to  this 
mathematical  model.  Also,  the  computer  solution  was  devised  in  such  a  manner  that  a 
wide  variety  of  flow  situations  can  be  considered;  this  depends  only  on  the  data  which  is 
provided  as  input  by  the  user.  However,  the  output  of  this  program  should  be  validated 
by  experimental  data  since  the  model  contains  many  assumptions. 
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Figure  6.  Calculated  Radial  Co-ordinate  of  Outside  Particle  as  a  Function  of  its 
Inertia  Parameter,  for  the  Jet  Width  Resulting  when  AF  =  1.0 


Calculated  Radial  Co-ordinate  of  Outside  Particle  as  a  Function 
tertia  Parameter,  for  the  Jet  Width  Resulting  when  AF  -  0.5 
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